Dynamics of the sub-Ohmic spin-boson model: a time-dependent variational study 
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The Dirac-Frenkel time-dependent variation is employed to probe the dynamics ol the zero tem- 
perature sub-Ohmic spin-boson model with strong friction utilizing the Davydov Di ansatz. It 
is shown that initial conditions of the phonon bath have considerable influence on the dynamics. 
Counterintuitively, even in the very strong coupling regime, quantum coherence features still man- 
age to survive under the polarized bath initial condition, while such features are absent under the 
factorized bath initial condition. In addition, a coherent-incoherent transition is found at a critical 
coupling strength a ~ 0.1 for s — 0.25 under the factorized bath initial condition. We quantify how 
faithfully our ansatz follows the Schrodinger equation, finding that the time-dependent variational 
approach is robust for strong dissipation and deep sub-Ohmic baths (s <^ 1). 
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I. INTRODUCTION 

It is of fundamental importance to study macroscopic 
behavior of open quantum systems under the influence 
of dissipative baths that they are inevitably in contact 
with [1, iH . As an open quantum system and its ther- 
mal bath with a finite number of modes together form 
an isolated system which conserves the energy, accord- 
ing to the Poincare recurrence theorem Q, the quan- 
tum system will eventually return to a state very close 
to its initial state. To circumvent this recurrence, the 
bath has to be expanded to contain an infinite number of 
phonon modes. The simplest model of a quantum system 
is a two-level system Hg — -~Aax/2, and Caldeira and 
LeggettQ have shown that a bath of harmonic oscillator, 
Hb = '^i^iblh, provides a very good approximation to 
real dissipative systems. Here ax is a Pauli matrix, and 
bj (bi) is the boson creation (annihilation) operator. In 
general, it is assumed that the two-level system is linearly 
coupled to the bath, Hsb = —o'z X^z ^li^l + In the 

absence of the bath, the quantum two-level system will 
oscillate between the states |-|-) and |— ) with frequency 
A (|±) are the two eigenstates of CTz), a quantum phe- 
nomenon which has no classical analog. Putting together 
all three Hamiltonian terms, 



H = -—ax 
2 



^c.,6t^,, + ^^A,(6j+60. (1) 
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The combined Hamiltonian described by Eq. ([T]) is called 
the spin-boson model. 

Despite its simplicity, the spin-boson model has been 
widely discussed in condensed phase physics and chem- 
istry, ranging from the process of electron transfer [5| to 
quantum entanglement ^ between a qubit and its 



bath. The coupling between the two-level system and 
the harmonic bath is completely specified by the spec- 
tral function 
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with a cutoff frequency ujc and a dimensionless constant 
a measuring the strength of the coupling. A bath de- 
scribed by a spectral function with s = 1 is referred to 
as an Ohmic bath. Both static and dynamical proper- 
ties of the spin-boson model with an Ohmic bath are 
well understood P, 0]. It is found that there exists a 
quantum phase transition from a non-degenerate delo- 
calized phase to a doubly degenerate localized phase and 
a turnover from a coherent phase to an incoherent one. 
The critical coupling for the quantum phase transition is 

« 1, and the critical coupling for the turnover to occur 
is aci ~ 0.5. In the Ohmic regime, the spin-boson model 
can be readily mapped onto the anisotropic Kondo model 
using bosonization techniques [l|, and known results can 
be borrowed from the Kondo model. However, success 
still eludes us in arriving at a correct description of the 
spin-boson physics in the sub-Ohmic regime (0 < s < 1). 

Several sophisticated numerical methods have been 
used to study the sub-Ohmic spin-boson model. An in- 
complete list includes: The numerical renormalization 
group method developed by Wilson ^ which reveals a 
continuous quantum phase transition for all < s < 1 
and weakly damped coherent oscillations on short time 
scales in the localized phase numerically exact real- 
time path integral method with quasiadiabatic propaga- 
tor revealing effective dynamic asymmetry in the pres- 
ence of a sub-Ohmic bath [l^, [ll|; quantum Monte 
Carlo method which determines the critical exponents 
for s < 1/2 [l3|; the sparse polynomial space representa- 
tion method which is based on the exact diagonalization 
to obtain numerical results [isj: the real-time path inte- 
gral Monte Carlo techniques which show that the coher- 
ent phase exists even in strong dissipation for s < 1/2 
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[13, [Hi; the numerically exact multilayer multiconfigu- 
ration time-dependent Hartree method (ML-MCTDH) 
which shows that the transition of the dynamics from 
weakly damped coherent motion to localization upon in- 
crease of the system-bath coupling strength [l^ [l3l ■ 

Unlike in the Ohmic regime, dynamics of sub-Ohmic 
spin-boson model is very sensitive to the initial condi- 
tions. There are two initial conditions of interest: one is 
the factorized initial condition with the bath in its vac- 
uum state initially; the other is the polarized initial con- 
dition consistent with typical experimental scenarios [2, 
under which the system is prepared in the ground state 
of Hb + HsB\a^=i- Many recent studies use the polarized 
initial condition [lTl4l3l [Tsj . The typical time scale of 
the spin dynamics is 1/A, the time it takes to hop from 
a spin state to another. On the other hand, the reorga- 
nization energy [isj . which describes the change in pop- 
ulation disparity between the two states |-t-) and |— ) as 
one goes from one initial condition to the other, is given 
as Jp dioJ{uj)/io = 2auicT{s), where r(s) is the gamma 
function of s. Ref. [l^ shows that the difference in spin 
dynamics under the two initial conditions is negligible 
for the Ohmic bath in the scaling limit A ^ Wc. How- 
ever, reorganization energies in the sub-Ohmic regime 
are larger than those in the Ohmic regime. It implies 
that much more time is needed for a sub-Ohmic bath 
to return to thermal equilibrium, and therefore the spin 
dynamics corresponding to the two initial conditions will 
display sizeable differences for certain parameter space 
(A/cjc,s). Ref. [m has confirmed such a physical pic- 
ture. It is commonly accepted that the increasing of the 
spin-bath coupling will eventually turn quantum coher- 
ent oscillations into classical-like damping, a picture sup- 
ported by the aforementioned numerical approaches for 
1/2 < s < 1. However, Ref. [l^ has recently claimed that 
such a picture may not always hold. Numerical data in 
the strong coupling regime, for example, with a = 30ac, 
show that the coherent phase exists for exponents up to 
s = 0.49. Such an "anomalous" result warrants further 
investigations. 

In this work, we adopt a variational framework to 
study the zero temperature dynamical properties of the 
sub-Ohmic spin-boson model. This is motivated by two 
facts: (i) The spin-bath interactions are formally iden- 
tical to the exciton-phonon coupling in a quasiparticle 
named a polarqn, which is generally described by the 
Holstein model [l9[ 

-^Holstcin ^ ^ ^Tn^ly^Q^m ^" ^ ^ Jmn{,^\',i^n ^ ^J^'^m) 

qm q 

where al^ (a,„) is the exciton creation (annihilation) op- 
erator, J,„„ is the hopping integral, and gq labels the 
exciton-phonon coupling strength; (ii) A time-dependent 
variational approach based on the Davydov ansatze has 
been widely used for describing the dynamics of Hol- 



stein systems. As a semi-classical approach for studying 
energy transport in deformable molecular chains, those 
ansatze (20| were put forward by Davydov and coworkers 
as two trial wave functions, namely, the Davydov Di and 
D2 ansatze. The first of Davydov's ansatze has the form 

|Di(i)) = Xa„(t)at |0)exexp[X(A„,(t)&t _ /i.e.)] |0)ph, 

71 q 

(4) 

where an(t) and A„g(t) are the variational parameters 
representing the exciton amplitude and the phonon dis- 
placements at the nth site, respectively, and |0)ox and 
|0)ph denote the exciton and the phonon vacuum states, 
respectively. The second of Davydov's ansatze is given 
by 

|D2(t)) = 5]«„(t)at |0)exexp[X(/3,(t)6t - /i.c.)]|0)ph, 

n q 

(5) 

Note that the phonon-displacement parameter in the 
much simplified \D2{t)), Pq{t), is independent of the 
site index n. The Davydov ansatze and their variants 
have also been applied successfully to study the one- 
dimensional Holstein polaron by Zhao and coworkers (2TI - 
[23 |. To probe polaron dynamics of the Holstein sys- 
tem using Davydov's ansatze [12, [13], we have employed 
the Dirac-Frenkel time-dependent variational principle, a 
powerful technique to obtain approximate dynamics for 
qimntum systems for which exact solutions are elusive 

By exploiting the analogy between the spin-boson 
model and the Holstein molecular crystal model, we pro- 
pose a time-dependent trial wave function very similar to 
the Davydov Di ansatz, and seek to develop an accurate 
description for dynamical properties of the sub-Ohmic 
spin-boson model with < s < 1/2 under both the po- 
larized and factorized initial conditions. It was pointed 
earlier [1^ that the Davydov ansatze bear close resem- 
blance to a multiconfigurational ansatz which contains 
more than one Slater determinant. We note that the 
ansatz employed in this work shares many characteris- 
tics with the multi-configurational Gaussian wave pack- 
ets (27l - [3l| , used in a variant of the powerful multicon- 
figuration time-dependent Hartree technique (MCTDH) 
[33 |. also known as the G-MCTDH method. Proposed in 
Ref. [13 and developed further in Refs. [li,[3l, the G- 
MCTDH method extends MCTDH to higher-dimensional 
systems by including a moving basis of Gaussian func- 
tions while restricted to an optimally chosen subset of 
the bath modes. It has been successfully applied to de- 
scribe dynamics of the Henon-Heiles potential model (ssj 
and the oscillator-bath model (28} . 

In the Holstein molecular crystal model, a similar reor- 
ganization energy [22j can be calculated from the phonon 
spectral density function using du}J{u})/uj, and is 
known to be proportional to the Huang- Rhys factor [13] , 
which measures the exciton-phonon coupling strength. 
For the factorized initial conditions, our previous stud- 
ies on the Holstein model reveal that the Davydov ansatz 
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arc especially accurate in the strong exciton-phonon cou- 
pling regime |2l|, a fact that will also be confirmed by 



where 



calculating the relative deviation for the sub-Ohmic spin- 
boson model studied here [22|. Furthermore, for a given 
spin-bath coupling strength, a smaller s yields a smaller 
relative deviation, inferring the highest accuracy of our 
ansatz in the deep sub-Ohmic regime s <C 1. A similar 
trend is believed for the polarized initial conditions. 

The paper is organized as follows. In Sec. II, we pro- 
pose a Davydov-like ansatz, and use the Dirac-Frenkel 
variational method to obtain the equations of motion for 
its parameters. In Sec. Ill, we present the numerical re- 
sults, which show the quantum coherence and entangle- 
ment for spin-boson model. The effect of different initial 
conditions is discussed in this section. In Sec. IV, con- 
clusions are drawn. 



II. METHODOLOGY 

Note that the spin-boson model can be viewed as a two- 
site Holstcin model with an infinite number of phonon 
modes in the one-exciton subspace. This equivalence be- 
tween the exciton-phonon coupling and the spin-boson 
interaction naturally leads to a trial wave function simi- 
lar to Davydov Di ansatz for the spin-boson model 

\D,{t)) = A(t)|+)exp[^(/Kt)&J - h.c.)]\0)ph 

+ B{t)\-)cxp[J2igi{t)bl ~ /i.c.)]|0)ph,(6) 



where A{t) and B{t) are complex variational parameters 
representing occupation amplitudes in states \+) and |— ), 
respectively, and //(t) and gi{t) label the corresponding 
complex phonon displacements of the 1-th phonon mode. 
In this work, we choose the Lagrangian formalism of the 
Dirac-Frenkel variational principle to obtain equations of 
motion for the variational parameters. The Lagrangian 
associated with the trial state \Ds{t)) is given as 



L^{Dsit)\'-—-H\Dsit)) 



(7) 



Substituting the trial state \Ds{t)) into the Lagrangian 
([7]). we arrive at the Lagrangian for the spin-boson model 



L= -{A*A- A*A + B*B - B*B) 



(^A*BgI2ifi'9l _|_ y^^*gE( /i9r)g-2 Eid/il +\9i\ ) 



'2 
1 

2 



lj2M\A\'ifi + fn-\B\'i9i+gt)] 



(9) 



The Dirac-Frenkel time-dependent variational principle 
yields the equations of motion for A, B, fi and gi 



d f dL\ dL 



dt V / 9u* 



0, 



(10) 



where w* denotes the complex conjugate of variational 
parameters m„, which can be A,B,fi or gi. From 
Eq. ([TU)) . one arrives at the equations of motion for A{t) 
and B{t) 

I I 

+ Mei:a/r9<-i(i/<i^+i«n] _ aEc.,|./zP , (11) 

= + *"f J2^9i9! - ft* 5/) + 9i) 



+ ^e^AaU.-U\h\'+\m\')] ^ BY,^i\gi\^ ■ (12) 



Similarly, the equations of motion for //(t) and 9i{t) are 
given as 

iAfi - - Auifi 
= ^(/,-.90e^a/r.-i(l/n=+l.<n] , (13) 
A; 

iBgi + B— - BuJigi 
^ M(gj_/;)e5:a9r/.-|(l/H^+l9<P)] . (14) 

Eqs. pT|) and ((T^ have been made use of to deduce the 
equations of motion for fi{t) and gi{t). It is found from 
the equations (fTTj) and (|12l) that 

±i\Af + \B\') = ^^{D,MDs{t))^0. (15) 

That the sum of \A\'^ and is conserved follows from 
early assignments of A{t) and B{t) in Eq. Therefore, 
+ which is the norm of \Ds{t)), can be set to 
unity 



\A\^ + = 1 



(16) 



{DS)\H\D,it)), 



(8) 



The equations of motion ((TTHH)) give a complete descrip- 
tion of the time evolution of \Ds{t)), and therefore, the 
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dynamics of the spin-boson model. In the spin-boson 
model, physical observables of interest are 

P,{t) = {a,) = {D,{t)\a,\Ds{t)) , i = x,y,z. (17) 

Here Px{t) describes the coherence between the |-t-) and 
— ) states, and Pz{t), the population difference. Upon 
substitution of the trial wave function ([6]) into Eq. ([TTl) , 
we obtain 

P,(t)= \A\^-\B\^ (18) 

_|_ ^5*gEi[9r/i-5(l/!l^ + l9il^)]^ (2.9) 

Py{t) = - i^*Be^'^'^""^(l-^'l'+l»'l')] 

+ iAB*e^'[3'^'-5(l-'''l'+lf'l')]- (20) 

Due to the invariance of H under ay — > —ay, one usu- 
ally has Py ~ for the ground state or thermal averages. 
However, the time-dependent observable Py(t) is in gen- 
eral nonzero. 

As mentioned in the Introduction, the initial condition 
has a vital influence on the dynamics of the spin-boson 
model with the sub-Ohmic bath. We assume that the 
spin is prepared in state |-|-) at t = 0, or A{0) = 1 and 
-8(0) = 0. At zero temperature, the factorized bath ini- 
tial condition corresponds to a phonon vacuum state with 
//(O) = gi{Q) = 0, while the polarized bath initial condi- 
tion is one in which the bath phonons arc in a displaccd- 
oscillator state with //(O) = .9/(0) = — A;/2a;;. 

When the phonon bath is absent, the dynamics of Pz (t) 
will be fully coherent and has no classical component. 
The bath will induce decoherence, and for a quantum 
dissipative system such as the spin-boson model, its pop- 
ulation difference has the form that Pz{t) ^ cos{^lt)e~'^^ 
[U, 0] on certain time scales. The oscillation cos(f]t) rep- 
resents the quantum coherence. The exponential decay 
is classical friction effect that induced by the bath. The 
dynamics is said to be coherent if $7 7^ 0, otherwise it is 
incoherent. On the other hand, the steady state is said to 
be localized if Pz{t 00) 7^ 0, otherwise it is delocalizcd 

Pz{t ^ 00) = o^M,M,Mm■ 

Another physical quantity of interest is the entangle- 
ment between the spin and the surrounding bath de- 
scribed by the von Neumann entropy S, also known as the 
entanglement entropy [13, ■ In the spin-boson model, 
it is given as @, [11] 

S = — w+ In w+ — w_ In , (21) 

where 

= {i±,/p[7p[+p!)/2 

= i±iy'l + 4|^P|B|2(e-E,l/<-.<P„l).(22) 

From Eq. (f20| . it is clear that Py{t) = if and only if 
A*Be^i ^' ^' is a real number, a condition that is satisfied 
in the ground state (or ensemble averages) due to Hamil- 
tonian invariance the transformation ay —ay m. 



III. NUMERICAL RESULTS 

The spectral function ([5]) gives full information for the 
spin-bath coupling A;. Together with equations of mo- 
tion pTI - fM)) . the dynamics of spin-boson model could 
be deduced for given specific initial conditions. The |-|-) 
state is usually populated at i = 0, i.e. A{0) = 1 and 
B{0) — 0. The initial conditions of the phonon bath are 
MO) = giiO) = and /i(0) = gi{0) = ~Xt/2oji for fac- 
torized and polarized initial conditions, respectively. We 
have to solve the 2Ni, + 2 equations of motion Eqs. (|11| - 
[T4| . together with the 2A^;, + 2 initial conditions men- 
tioned above, where Nh is the number of the phonon 
modes considered. 
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FIG. 1: Under the polarized bath initial conditions, the time 
dependent population difference Pz{t) for s = 0.25 (ac ~ 
0.022) and A/cjc = 0.1 is presented. Seven values of a are 
taken. The inset is a magnified figure for large couplings. 

We will adopt the homogeneous discretization proce- 
dure used in Ref. [i^. The frequencies of the Nb har- 
monic modes are equally distributed in the frequency 
range w G [uj min, uj mux] with spacing Acj = ujmax/Nb 
so that ojmin = cji = Aw and ui ~ ll^u. The fre- 
quency spacing Aw determines the Poincare recurrence 
time Tp = 27r/Aw which must be greater than any time 
scale of interest [4l[, [4^ . Throughout this work, we will 
use Nb = 20000 and Wmax = 4a;c, resulting in a recur- 
rence time Tp = lOOOOTr that places our simulations at a 
safe distance from the Poincare recurrence. Correspond- 
ingly, from the integration of the spectral density over 
w 

Nb poo Nb 

^Af=/ dujj{uj) ^^J{uji)Auj, (23) 
1=1 /=i 

we obtain that A^ = J{uJi)Au!. It is found that un- 
der the factorized bath initial condition the simulation 
results are insensitive to the number of phonon modes. 
Under the polarized initial condition, the number of the 
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deviation [22] of the trial state \Ds{t)) defined as 



FIG. 2: Time dependent population difference Pzit) (upper 
panel) and relative deviation a{t) (lower panel) for s = 0.25 
(ac ~ 0.022) and A/tOc = 0.1 under the factorized bath initial 
conditions. Seven values of a are used for comparison. 



phonon modes has a considerable influenee on the dy- 
namics, which is especially true in the strong coupling 
regime. However, our numerical tests show that good 
convergence is reached when Nt = 20000 for the time 
periods considered. 

As shown in Fig. 1, under the polarized initial con- 
dition, the population difference as a function of time, 
Pzit), manifests coherent oscillations even for very large 
couplings a = 0.3 sa ISctcj where ac ~ 0.022 [TTl |43| is 
the critical coupling for the quantum phase transition, 
ft is clearly seen that oscillatory behavior emerges even 
for very strong coupling far beyond ac- Furthermore, 
the characteristic oscillation frequency of Pz [t) increases 
with increasing a. Our results agree with those of Kast 
et al. obtained u sing the real-time path integral Monte 
Carlo simulation [l^. It is widely accepted that under 
the polarized initial condition, a quantum dissipative sys- 
tem such as the spin-boson model is expected to display 
classical over-damped behavior (or incoherent phase) at 
strong spin-bath couplings. Our results reveal that this 
is not the case for < s < 0.5. However, such apparent 
contradictions only appear for the polarized bath initial 
condition. For the factorized initial bath condition, the 
persistent coherence does not occur. The upper panel 
of Fig. 2 shows the population difference as a function 
of time, Pz{t), under the factorized bath initial condi- 
tion. It is found that the critical coupling strength for 



the coherent-incoherent transition is a 



(/) 

CI 



0.1 « A.bar 



a{t) 



vmm) 



bath 



(24) 



where i?bath denotes the average energy of bath within 
the time interval considered, and \5{t)) is the devia- 
tion vector quantifying how faithfully \Ds{t)) follows the 
Schrodinger equation: 



\m) = {^dt-H)\Ds{t)). 



(25) 



In another word, the smaller the relative deviation (T(t), 
the closer the trial state \Ds{t)) obeys the Schrodinger 
equation. To compare with the ML-MCTDH method, a 
calculation is carried out for four values of s. The up- 
per panel of Fig. 3 shows the time-dependent population 
difference Pz{t) for a = 0.2, Wc/A = 5, in good agree- 
ment with corresponding results in Ref. (l7j . As shown 
in the lower panel of Fig. 3, the relative deviation a{t) 
gradually increases with time before reaching a satura- 
tion value that have a strong dependence on the exponent 
s, indicating that the smaller the exponent s is, the more 
accurate our ansatz becomes. 

Thus, the bath initial conditions play an important 
role in the dynamics of the sub-Ohmic spin-boson model. 
In Ref. [isj . the authors studied out-of-equilibrium bath 
preparations with respect to the initial state of the spin. 
The polarized bath initial condition corresponds to the 
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In the lower panel of Fig. 2, we also plot the relative 



FIG. 3: Time-dependent population difference P2(t) (upper 
panel) and relative deviation a(t) (lower panel) for a = 
0.5, uJc//S. = 5, and four values of the exponent s: s = 
0.25, 0.5, 0.75, and 1. The factorized bath initial condition 
is taken. 
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FIG. 4: Population difference Pz{t) as a function of time for 
s = 0.25 (oc ~ 0.022) and A/ljc = 0.1 with polarized and 
factorized bath initial conditions. 



case where the bath distribution is fully equilibrated with 
the initial state of the spin, while the factorized bath ini- 
tial condition corresponds to the one most displaced from 
the equilibrium. We note that our results for a = 0.1 in 
Fig. 1 and in the upper panel of Fig. 2 are consistent with 
those in Ref. [l^, where it is found that compared with 
the polarized initial condition, the factorized one yields 
decreased oscillation frequency of the dynamics and in- 
creased initial loss in population (cf. the pink lines in 
Fig. 1 and the upper panel of Fig. 2). 

An explicit comparison between the polarized and fac- 
torized initial conditions is given in Fig. 4, where sub- 
stantial difFerenccs are revealed between the two bath 
initial conditions for a = 0.9q;c and a ~ 1.3ac. Un- 
der the factorized bath initial condition, the oscillations 
occur around zero average in both the localized and delo- 
calized phases, and the population difference Pzit) is in 
a delocalized phase even for a > ac, a parameter regime 
where the state is expected to be localized. In contrast, 
the oscillations occur around finite values even in the de- 
localized phase for the polarized initial condition. Also, 
the damped constant (or steady state) is an increasing 
function of a. Similar results are also found by Nalbach 
and Thorwart using the quasiadiabatic propagator path 
integral. [lT| . 

We also monitor the entanglement between the spin 
and the bath via the von Neumann entropy. At i = 0, the 
systems is in a separable state, so that 5'(0) = 0. Eq. ([2T|) 
manifests that the entropy S increases monotonously 
with w-f from and reaches its maximum at a;+ = 1/2, 
then decreases monotonously to 0. Figs. 5 and 6 show 
that there are substantial differences in the time evo- 
lution of the entropy under the polarized and factorized 
bath initial conditions. Under the polarized initial condi- 
tion, overall the entropy decreases as a increases as shown 
in Fig. 5. For strong coupling strengths, the entropy 
eventually vanishes as expected Q. Under the factorized 
initial condition, the entropy establishes its steady val- 



ues quickly for various coupling strengths. Interestingly, 
the steady value is not a monotonous function of a, and 
reaches its maximum at approximately a w 0.07. 



IV. CONCLUSIONS AND DISCUSSIONS 

Although quite a few numerical approaches have been 
applied to study the ground state and dynamical prop- 
erties of the sub-Ohmic spin-boson model, few analyt- 
ical treatments have been available. Recently, Chin et 
al. [1^ used an extension of the Silbey- Harris variational 
wave function to study the ground state of the sub- 
Ohmic Spin-boson model with < s < 0.5, and found 
that such a trial state generates correct mean-field expo- 
nents for the continuous localization-delocalization tran- 
sition. The asymmetrically displaced-oscillator (ADO) 
trial state used in Ref. [i^ is of the form 

|^)ado = A|+)exp[^(/i6t _ /i.c.)]|0)ph 

I 

+ i?|-)exp[^(.9(&t_/,.c.)]|0)pi„ (26) 
I 

where the variational parameters A, B, f^s and g'^s are 
all real numbers to be determined by the ground-state 
energy minimization. It is interesting to note that our 
ansatz, Eq. is reduced to Eq. ([26)) if we restrict all 
variational paremeters to be real, and consequently, the 
Dirac-Frenkel variation employed in this work is reduced 
to the conventional variational principle for the ground 
state used in Ref. [i^. Therefore, we have successfully 
extended the static trial state of Eq. ([26| to its dynamical 
counterpart, which is similar in form to the Davydov Di 
ansatz. The foregoing connection also helps map out 
the validity regime of our ansatz. As pointed out by 
Nazir et al, although the ADO state works well for the 
sub-Ohmic baths with s < 0.5, it becomes unstable and 




FIG. 5: Von Neumann entropy calculated with the polarized 
bath initial condition for s — 0.25 (oc ~ 0.022) and A/cjc = 
0.1. 
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FIG. 6: Von Neumann entropy calculated under with the fac- 
torized bath initial condition for s — 0.25 (qc ~ 0.022) and 
A/lJc = 0.1. 

deviates from the well-established results for the Ohmic 
case and stron g co uphnR (the cause of which is still under 
investigation) |44| . Thus, it is expected that our ansatz 
reveals reliable dynamics of the sub-Ohmic spin-boson 
model in the regime of s -C 1, but may lose accuracy for 
large s. 

For the sub-Ohmic spin-boson model with < s < 1/2, 
detailed dynamics of the spin-boson model in the strong 
coupling regime is still surrounded by contention. The 
hierarchy of the Davydov ansatze of varying sophistica- 
tion has been known to be competent in handling polaron 
dynamics in the strong coupling regime [l^, [l^l • Using a 
version of the most accurate of the hierarchy, the Davy- 
dov Di ansatz, we have carried out a time-dependent 
variational calculation with regard to the dynamic prop- 
erties of the sub-Ohmic spin-boson model. It is a sim- 



ple, yet extremely efficient approach to investigate the 
dynamics of a quantum dissipativc system, such as the 
population disparity Pz{t) of the sub-Ohmic spin-boson 
model. Our results are consistent with those obtained 
using numerically much more expensive advanced nu- 
merical methods, such as the path integral Monte Carlo 
simulations [III [lal and the ML-MCTDH technique [l7j . 
It is found that the bath initial conditions have consider- 
able influence over the dynamics of this many-body dissi- 
pativc system. Even in the very strong coupling regime, 
quantum coherence features still manage to survive under 
the polarized bath initial condition, while such features 
are absent under the factorized bath initial condition. 
The onset of the incoherent phase occurs at a = 0.1 for 
s = 0.25 under the factorized bath initial condition. Our 
findings arc consistent with those in Ref. [isj , which first 
reported the persistence of coherent quantum dynamics 
at strong dissipation under the polarized bath initial con- 
dition. Furthermore, the Davydov Di ansatz has been 
employed successfully to study excitation energy trans- 
fer across light-harvesting complexes in photosynthesis 
[23 |. Our approach may turn out to be a competitive 
tool to investigate sustained quantum coherence recently 
discovered to reside in pigment networks even at elevated 
temperatures (4^-113] ■ 
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